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Abstract 

One-  and  two-parameter  numerical  optimization  analyses  of  the  cathode  catalyst  layer  of  the  PEM  fuel  cell  were  investigated  with  the 
objective  of  optimizing  the  current  density  of  the  catalyst  layer  at  a  given  electrode  potential.  Catalyst  design  parameters,  such  as  Nafion 
content,  platinum  loading,  catalyst  layer  thickness  and  porosity,  were  considered.  Numerical  analysis  shows  that  there  is  a  global,  optimal 
solution  for  each  one-parameter  optimization.  A  global,  optimal  solution  also  exists  for  two-parameter  optimizations.  The  optimal  solution 
and  the  surface  of  the  local  current  density  as  a  function  of  two  design  parameters  are  given  for  each  two-parameter  optimization.  The 
numerical  analyses  provided  here  are  useful  in  the  design  of  cathode  catalyst  layers. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  are  promis¬ 
ing  candidates  for  next  generation  power  sources.  They 
have  received  much  attention  in  recent  years  due  to  their 
“zero-emission”  quality  and  relative  high  power  density.  A 
PEM  fuel  cell  generally  consists  of  a  membrane  electrode 
assembly  (MEA),  sandwiched  between  two  bipolar  plates 
that  allow  the  introduction  of  separate  reactant  gases.  The 
MEA  consists  of  a  proton  conducting  electrolyte  mem¬ 
brane  sandwiched  between  two  porous  gas  diffusion  layers 
(GDLs)  which  have  been  coated  on  one  side  with  a  thin 
catalyst  layer:  it  is  often  prepared  from  a  mixture  of  catalyst 
dispersed  on  carbon  black,  polyelectrolyte,  and  a  hydropho¬ 
bic  agent.  Hydrogen  gas  diffuses  through  the  anode  GDL  to 
the  anode  catalyst  layer  and  is  oxidized  to  protons.  Concur¬ 
rently,  the  electrons  travel  through  the  electrical  circuit  to  the 
cathode  where,  in  the  presence  of  protons  that  have  diffused 
through  the  membrane,  they  reduce  oxygen  that  has  diffused 
through  the  cathode  GDL  to  the  cathode  catalyst  layer. 

Many  different  transport  processes  occur  in  different  re¬ 
gions  of  the  PEM  fuel  cell.  These  include  mass  transfer  of 
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the  chemical  species  (typically  hydrogen  and  oxygen)  in  the 
gas  channels  and  porous  GDLs,  proton  and  electron  trans¬ 
port  through  the  catalyst  electrodes  and  membrane,  water 
transport  in  different  phases  through  catalyst  layers,  GDLs 
and  the  membrane,  and  heat  transfer.  Among  these  different 
components  and  different  transport  processes,  the  transport 
processes  in  the  cathode  catalyst  layer  play  a  significant  role 
in  the  performance  of  the  PEM  fuel  cell  because  this  is  the 
region  where  the  rate  limiting  chemical  reaction  takes  place 
and  potential  power  is  wasted.  Therefore,  research  on  cath¬ 
ode  performance  has  been  very  active,  and  many  models 
have  been  proposed  [1-5]. 

Commercialization  requires  that  the  fuel  cell  achieves  its 
best  performance  at  the  lowest  cost.  Optimal  PEM  fuel  cell 
performance  should  be:  the  highest  power  output,  highest 
fuel  efficiency,  while  maintaining  reliability  and  durability, 
and  at  the  lowest  cost.  Assuming  that  adequate  water  and 
heat  management  can  be  achieved,  optimal  performance  will 
strongly  correlate  with  an  optimized  cathode  catalyst  layer. 

Optimization  by  experiment  [6-8]  is  time  consuming  and 
costly  compared  to  numerical  optimization  [9,10],  in  partic¬ 
ular,  when  more  than  one  variable  is  involved.  On  the  other 
hand,  the  numerical  optimization  method  is  inexpensive  and 
can  yield  global  optimization  because  many  parameters  can 
be  systematically  varied.  Qi  and  Kaufman  [8]  investigated 
experimentally  the  effects  of  Nafion  content  and  platinum 
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specific  reaction  surface  area  per  volume  of  the  catalyst  layer  (cm-1) 
catalyst  surface  area  per  unit  mass  of  the  catalyst  (cmg-1) 
oxygen  concentration  in  the  catalyst  layer  (mol  cm-3) 
reference  oxygen  concentration  (mol  cm-3) 

oxygen  concentration  at  the  interface  of  the  GDL  and  the  catalyst  layer  (mol  cm-3) 

bulk  diffusion  coefficient  of  oxygen  in  electrolyte  (cm2  s-1) 

bulk  diffusion  coefficient  of  oxygen  gas  (cm2  s-1) 

diffusion  coefficient  of  oxygen  in  Nation  (cm2  s-1) 

diffusion  coefficient  of  oxygen  in  liquid  water  (cm2  s-1) 

effective  diffusion  coefficient  for  oxygen  through  the  catalyst  layer  (cm2  s-1) 

effective  diffusion  coefficient  of  oxygen  in  GDL  (cm2  s-1) 

open  circuit  voltage  (V) 

Faraday’s  constant  (Cmol-1) 

reference  exchange  current  density  (A  cm-2) 

proton  current  density  at  x  (A  cm-2) 

current  density  at  a  given  electrode  potential  Vo  (A  cm-2) 

total  current  density  in  the  catalyst  layer  (A  cm-2) 

Henry’s  law  coefficient  (dimensionless) 
thickness  of  the  cathode  catalyst  layer  (cm) 
thickness  of  the  GDL  (cm) 

catalyst  mass  loading  per  unit  area  of  the  catalyst  layer  (gem-2) 

number  of  electrons  transferred  in  cathode  reaction 

oxygen  pressure  (atm) 

gas  constant  (JK-1  mol-1) 

temperature  (K) 

given  electrode  potential  (V) 

coordinate  (cm) 

mass  percentage  of  Nation  to  the  sum  of  Nation  and  solid  particles 
mass  percentage  of  platinum  catalyst  on  the  support  carbon  black 
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cathodic  transfer  coefficient 
anodic  transfer  coefficient 

volume  fraction  of  the  Nation  in  the  catalyst  layer 
volume  fraction  of  solid  catalyst  particles  in  the  catalyst  layer 
volume  fraction  of  the  void  space  in  the  catalyst  layer, 
porosity  of  the  GDL 
overpotential  (V) 

bulk  conductivity  of  Nation  (S  cm-1) 

bulk  conductivity  of  solid  catalyst  (platinum  +  carbon)  (Scm-1) 

mass  density  of  the  platinum  (gem-3) 

mass  density  of  the  carbon  black  (gem-3) 

mass  density  of  Nation  (gem-3) 

effective  conductivity  in  Nation  (S  cm-1) 

effective  conductivity  in  solid  catalyst  particle  (S  cm-1) 


loading  on  cathode  performance  individually  and  concluded 
that  30wt.%  of  Nation  content  and  0.20  db  0.05  mg/cm2  of 
platinum  loading  were  the  optimum  values  to  obtain  the  best 
performance  when  20%  Pt/C  catalyst  was  used  to  make  the 
cathode.  The  reported  best  performance  of  30wt.%  Nation 
or  0.20  ±  0.05  mg/cm2  Pt  loading  seems  to  be  local  maxi¬ 


mum,  not  global.  They  did  not  discuss  what  values  the  opti¬ 
mal  Nation  content  and  platinum  loading  would  be  in  order 
to  obtain  the  best  cathode  performance  when  both  Nation 
content  and  platinum  loading  are  variables.  Boyer  et  al.  [9] 
gave  a  set  of  equations  to  optimize  the  thickness,  catalyst 
loading  and  Nation  content  individually  when  only  one  mass 
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transfer  process  was  significant,  and  the  numerical  analysis 
and  the  experimental  results  were  compared.  Marr  and  Li 
[10]  performed  a  numerical  analysis  to  optimize  the  perfor¬ 
mance  of  cathode  and  indicated  that  optimal  ionomer  con¬ 
tent  as  well  as  an  optimal  void  fraction  exists  regarding  the 
change  of  the  local  current  density  with  ionomer  content  or 
void  fraction.  In  their  calculation,  the  effective  diffusion  co¬ 
efficient  of  oxygen  in  catalyst  layer  is  estimated  by  applying 
the  Bruggeman  relation  to  the  combination  of  ionomer  and 
water  in  which  the  bulk  diffusion  coefficient  is  calculated 
based  on  a  series  resistance  model. 

This  work  focuses  on  the  numerical  optimization  of  the 
cathode  catalyst  layer  with  respect  to  one  or  two  of  the 
four  design  parameters:  Nation  content,  void  volume  frac¬ 
tion  or  porosity,  thickness  and  platinum  loading,  by  means  of 
the  one-dimensional  “macro-homogeneous”  model  [3,10]. 
Isothermal  and  steady  state  assumptions  are  assumed.  Oxy¬ 
gen  diffusion  is  described  by  Fick’s  law.  The  effective  dif¬ 
fusion  coefficient  of  oxygen  in  catalyst  layer  is  estimated 
as  follows:  first  applying  the  Bruggeman  relation  to  Nation 
and  water  to  obtain  the  effective  diffusion  coefficients  of 
oxygen  in  Nation  and  water,  respectively,  and  then  using  a 
series  resistance  model.  Proton  migration  in  catalyst  layer 
is  described  by  the  Ohm’ s  law,  and  electrochemical  kinetics 
are  described  by  the  Butler- Volmer  equation. 


2.  Mathematical  model  for  the  catalyst  layer 


catalyst  layer,  n  the  number  of  electrons  transferred  in  the 
oxygen  reduction  reaction,  and  F  is  the  Faraday’s  constant. 
Dff  ,  the  effective  diffusion  coefficient  for  oxygen  in  the  cat¬ 
alyst  layer  is  evaluated  by  applying  series  resistance  model 
[10]  to  effective  diffusion  of  oxygen  in  Nation  and  liquid 
water: 
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where  N  and  w  are  the  effective  diffusion  coeffi¬ 
cients  of  oxygen  in  Nation  and  liquid  water,  respectively, 
and  £n  and  sy  are  the  volume  fractions  of  Nation  and  void 
space  in  the  catalyst  layer,  respectively. 

According  to  the  model  assumption  that  the  catalyst  layer 
is  fully  flooded  and  the  void  space  in  the  catalyst  layer  is 
completely  filled  with  water,  the  effective  diffusion  coeffi¬ 
cients  of  oxygen  in  Nation  and  water  can  be  evaluated  re¬ 
spectively  using  Bruggeman  relation: 


T'jeff  r\  3/2 

^o2,n  —  ^o2,N% 

(3a) 

and 

T^eff  r\  3/2 

^02,W  —  ^02,W£y 

(3b) 

where  £>o2,n  and  £>o2,w  are  the  bulk  diffusion  coefficients 
of  oxygen  in  Nation  and  water,  respectively. 

The  oxygen  concentration  in  the  catalyst  layer,  c o2,  is 
restricted  by  following  boundary  condition: 

co2(x  =  0)  =  c*  (4) 


A  schematic  diagram  of  the  cathode  catalyst  layer  is 
shown  in  Fig.  1.  Oxygen  diffusion  through  the  catalyst  layer 
is  assumed  to  be  due  to  its  concentration  gradient,  and  there¬ 
fore  the  oxygen  concentration  distribution  is  governed  by 
Fick’s  law: 

dcp2  _  /(*)  -  hot 
dx  ~  nFDf 

where  co2  is  the  concentration  of  oxygen  in  the  catalyst 
layer,  v  the  coordinate,  I(x)  the  proton  current  density  in  the 
catalyst  layer,  /tot  =  I(L)  the  total  current  density  in  the 


where  c*  is  the  oxygen  concentration  at  the  GDL/catalyst 
layer  interface.  Applying  Fick’s  law  for  oxygen  diffusion 
in  the  GDL,  and  solving  the  resulting  equation,  yield  the 
following  expression  for  c*: 

*  1  /  4otTGDL 

Ko2  y  RT  /iFD^  gdl 

where  Po2  is  the  oxygen  pressure,  R  the  gas  constant,  T 
the  temperature,  Lqdl  the  thickness  of  the  GDL,  and  Kq2 
is  Henry’s  coefficient  (dimensionless)  for  oxygen  gas  dis¬ 
solved  in  liquid  water  and  can  be  evaluated  from  the  empir¬ 
ical  relation  [3]  as 

1/666  \ 

Ko 2  =  —e*P  ( — jr  +  14.1 )  (6) 


Dq  gdl  is  the  effective  diffusion  coefficient  of  oxygen  in 
GDL  and  can  also  be  calculated  by  Bruggeman  relation: 

^o2,gdl  —  Ad2,gdl£Gdl  ( 7 ) 


where  Do2,gdl  is  the  bulk  diffusion  coefficient  of  oxygen 
gas  and  £qdl  is  the  porosity  of  the  GDL. 

Applying  Ohm’ s  law  to  both  the  proton  and  electron  trans¬ 
port  yields  the  governing  differential  equation  for  the  over¬ 
potential,  77,  in  the  cathode  catalyst  layer: 


drj  _  (  i  (  1  \  ^  1 

d~X  ~  +  ^eff  J  7«  -  ^eff  7‘°‘ 
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where  a^IT  and  a^IT  are  the  effective  proton  conductivity 
of  Nation  and  electronic  conductivity  of  the  solid  catalyst 
particles  (carbon-supported  Pt),  respectively,  and  estimated 
from  the  Burggeman  relation  to  account  for  the  porous  nature 
of  the  cathode  catalyst  layer: 


.eff 


3/2 


eff  _  ^  3/2 

aS  ~  as£s 


(9) 

(10) 


where  a n  and  as  are  the  bulk  proton  and  electronic  conduc¬ 
tivities  of  Nation  and  solid  catalyst  particles,  respectively. 
£s  is  the  volume  fraction  of  the  solid  catalyst  particles  in  the 
catalyst  layer  and  is  given  by  [10]: 
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where  ppt  and  pc  are  the  mass  densities  of  the  Pt  and  carbon 
black,  respectively;  %Pt  represents  the  mass  percentage  of 
Pt  catalyst  supported  on  the  carbon  black;  rapt  is  the  Pt  mass 
loading  per  unit  area  of  the  catalyst  layer. 

The  sum  of  the  Nation  volume  fraction  (£n),  the  void 
space  volume  fraction  (£y)  and  the  solid  catalyst  particle 
volume  fraction  (£s)  is  equal  to  1: 

£N  H"  £V  T  £S  =  1  (12) 


Combining  the  proton  mass  balance  and  the  Butler- Volmer 
equation  gives  the  equation  for  the  proton  current  density: 
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where  A  is  the  specific  reaction  surface  area  per  unit  vol¬ 
ume  of  the  catalyst  layer,  /gef  the  reference  exchange  cur¬ 
rent  density,  c{^  the  reference  oxygen  concentration,  ac  and 
a a  are  the  cathodic  and  anodic  transfer  coefficients,  respec¬ 
tively.  The  current  density,  I(x ),  should  satisfy  the  following 
boundary  conditions: 


7(0)  =  0;  7(L)  =  7tot 


(14) 


According  to  experimental  data  [11],  the  reference  exchange 
current  density  for  oxygen  reduction  in  Nafion  1 17  is  related 
with  the  cell  temperature  by  [3] 


j-ref  _  ^q(3.507-4001/7) 


(15) 


The  specific  reaction  surface  area  of  the  catalyst  layer  is  re¬ 
lated  to  the  Pt  loading,  rapt,  and  the  catalyst  layer  thickness, 
L,  as  following: 
rapt 

A  =  A0-^  (16) 

where  Aq  is  the  catalyst  surface  area  per  unit  mass  of  catalyst. 


3.  Statement  of  the  optimization  problem 

In  this  optimization  problem,  we  consider  the  following 
four  design  variables:  Nafion  volume  fraction  (£n),  void 


space  volume  fraction  (£y),  Pt  loading  (rapt)  and  catalyst 
layer  thickness  (L).  Since  PEM  fuel  cells  are  usually  oper¬ 
ated  in  a  specific  voltage  range,  the  objective  of  this  opti¬ 
mization  problem  is  to  maximize  the  current  density  at  a 
given  voltage. 

Denoting  the  current  density  generated  in  the  catalyst 
layer  as  7q  at  a  given  electrode  potential  Vo,  the  general  op¬ 


timization  problem  can  be  stated  as 

Maximize  7o  (17a) 

with  respect  to  £n,  £y,  L  and  rapt,  and  subject  to 
0  <  £n  <  1  (17b) 

0  <  £w  <  1  (17c) 

0  <  £S  <  1  (17d) 

and 

£n  +  £w  +  £s  =  1  (17e) 


where  £s  is  given  by  Eq.  (11). 

The  weight  percentage  of  Nafion  to  the  sum  of  Nafion  and 
solid  particles,  denoted  by  %N,  is  usually  used  to  measure 
the  Nafion  content  in  the  electrode  and  can  be  calculated 
from  the  volume  fraction  as 

PN^N 


%N  = 


Pn£n  +  (l/%Pt)  (rapt/ L) 


(18) 


where  is  the  mass  density  of  Nafion  in  liquid  water. 

When  one  parameter  varies  and  the  rest  three  are  fixed,  we 
have  four  one-parameter  optimizations  corresponding  to  the 
four  parameters:  £n,  £v,  L  and  rapt,  respectively.  When  two 
of  the  four  parameters  vary  and  the  other  two  are  fixed,  there 
are  six  two-parameter  optimizations.  In  this  paper,  we  will 
discuss  the  one-  and  two-parameter  optimization  problems 
to  optimize  the  catalyst  layer  performance. 


4.  Numerical  results  and  discussions 

Fig.  2  gives  the  comparison  of  the  experimental  data  [8] 
and  the  fitting  current  density  at  electrode  potential  Vo  = 
0.6  V.  The  optimization  shows  that  when  the  Nafion  con¬ 
tent  is  0.32087  (here  the  other  design  parameters  we:  L  = 
11.8  pan,  £y  =  0.10333,  rapt  =  0.332  mg/cm2),  the  current 
density  at  electrode  potential  Vo  =  0.6  V  reaches  its  maxi¬ 
mum  or  0.93583  A/cm2. 

Fig.  3  plots  the  experimental  data  [8]  and  the  fitting  cur¬ 
rent  density  curve  as  a  function  of  platinum  loading.  The 
highest  current  density  is  0.7 1492  A/cm 2  when  the  platinum 
loading  is  0.21032  mg/cm2  which  is  very  close  to  Qi  and 
Kaufman’s  conclusion  [8]:  an  optimum  platinum  loading  is 
around  0.20  db  0.05  mg/cm2. 

The  following  section  will  discuss  the  numerical  opti¬ 
mization  analysis  corresponding  to  one-  or  two-parameter 
optimization  problems.  The  parameter  values  (baseline  data) 
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Fig.  2.  Comparison  of  experimental  data  [8]  and  fitting  result  (Vo  =  0.6  V). 


used  in  analysis  are  given  in  the  Table  1  and  the  polariza¬ 
tion  curve  computed  using  these  data  is  shown  in  Fig.  4.  In 
the  process  of  optimization,  the  corresponding  optimization 
parameter  is  no  longer  regarded  as  a  fixed  value  but  as  a 
variable,  satisfying  the  corresponding  constraint  conditions. 

Fig.  5  gives  the  current  density  at  a  given  electrode  poten¬ 
tial  (Vo  =  0.6  V)  as  a  function  of  Nafion  volume  fraction. 
From  this  figure  it  can  be  easily  observed  how  the  Nafion 
volume  fraction  affects  the  current  density.  Fig.  5  indicates 
the  optimal  Nafion  volume  fraction,  £n,  is  0.33619  (here 
L  =  11.8  pan,  £n  =  0.08801  and  mpt  =  0.332mg/cm2), 
and  provides  a  maximum  current  density  of  0.37467  A/cm2. 
However,  the  asymmetric  broadness  of  the  plot  indicates 
that  although  the  current  density  is  not  very  sensitive  to  the 
volume  fraction  near  this  optimal  value  (from  0.3  to  0.35), 
it  is  more  sensitive  to  larger  values  of  Nafion  content  than 
to  lower  values  of  Nafion  content.  For  example,  the  current 


Fig.  3.  Comparison  of  experimental  data  [8]  and  fitting  result  for  Pt 
loading  (Vo  =  0.6  V). 


Table  1 

Parameter  values  used  in  the  baseline  calculation 


F  =  96485  (Cmol-1) 

R  =  8.315  (JK-1  mol-1) 
n  =  4 
ol  c  =  1 
aa  —  0.5 
E  =  1.23  (V) 

Poo  =0.21  (atm) 

T  =  308  (K) 

£>o2,gdl  =  2.396  x  10“ 1  (cm2  s-1) 
Z)o2,n  =  1-844  x  10“6  (cm2s_1) 
P>o2,w  =  3.032  x  10-5  (cm2  s_1) 
%Pt  =  0.2 

rapt  =  0.332  x  10“ 1  (gem-2) 


A0  =  112  x  104  (cmg-1) 

=  1.2  x  10-6  (mol  cm-3) 
L  =  1.18  x  10-3  (cm) 

Lgdl  =  100.0  x  10-2  (cm) 
V0  =  0.6  (V) 

£n  —  0.3 

£gdl  =0.3 

on  —  0.17  (S  cm-1) 

as  =  7.27  x  102  (Scm”1) 

pPt  =  21.5  (gem-3) 

pc  -  2.0  (gem-3) 

pN  -  2.0  (g  cm-3) 


data. 


Fig.  5.  Optimal  point  and  current  curve  as  a  function  of  Nafion  volume 
fraction  £n  (Vq  =  0.6  V). 
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Fig.  6.  Optimal  point  and  current  curve  as  a  function  of  void  volume 
fraction  s\  (V0  =  0.6  V). 


density  drops  by  10%  when  the  Nation  volume  fraction  falls 
below  0.23,  but  the  catalyst  layer  only  tolerate  up  to  a  vol¬ 
ume  fraction  of  0.4  before  the  current  density  falls  below 
10%  of  the  maximum. 

Fig.  6  shows  the  change  of  current  density  with  the  void 
volume  fraction  of  the  catalyst  layer.  The  optimal  void  vol¬ 
ume  fraction  £y,  is  0.087981  (here  L  =  11.8  pan,  £n  = 
0.33622  and  mpt  =  0.332  mg/cm2),  and  provides  a  current 
density  of  0.37467  A/cm2.  In  contrast  to  the  case  of  £ n,  the 
dependence  of  the  current  density  on  £y,  is  much  more  sensi¬ 
tive  to  lower  void  volume  fractions  than  to  higher.  The  max¬ 
imum  current  density  in  this  case  is  exactly  same  with  the 
value  of  the  case  £n  because  in  both  cases,  the  volume  frac¬ 
tion  of  solid  catalyst  particles  is  a  constant  which  means  the 
sum  of  the  void  volume  fraction  and  Nation  volume  fraction 
does  not  change  (1  —  £s),  so  the  current  density  as  a  func¬ 
tion  of  void  volume  fraction  can  be  regarded  as  the  result  of 
substitution  of  Nation  volume  fraction  by  (1  —  £s  —  £y)  in 
the  current  density  function  of  Nation  volume  fraction. 

Fig.  7  shows  the  relationship  between  current  density  and 
Pt  loading.  The  optimal  loading  is  0.24365  mg/cm2  (here 
L  =  11.8  jam,  £n  =  0.3  and  £y  =  0.27743),  and  provides 
a  maximum  current  density  of  0.39768  A/cm2.  The  current 
density  is  not  too  sensitive  to  the  Pt  loading  near  the  optimal 
value.  A  variation  of  about  ±45%  in  Pt  loading,  leads  to 
about  10%  decrease  in  calculated  maximum  current. 

The  calculated  current  density  is  even  more  sensitive  to 
thinner  thickness  of  the  catalyst  layer  than  to  thicker  one 
(Fig.  8).  Having  a  maximum  value  of  0.38606  A/cm2  at  op¬ 
timal  thickness  13.0972  pan,  a  variation  of  only  15%  less 
than  the  optimal  thickness,  leads  to  10%  decrease  in  calcu¬ 
lated  current  while  a  variation  of  60%  more  than  the  opti¬ 
mal  thickness  leads  to  10%  decrease.  Hence,  in  the  design 
of  catalyst  layers,  thickness  needs  special  consideration. 

Fig.  9  shows  the  optimal  volume  fractions  for  Nation 
and  void  space  in  the  two-parameter  optimization  when  the 


Fig.  7.  Optimal  point  and  current  curve  as  a  function  of  Pt  loading  mpt 
(Vo  =  0.6  V). 


Fig.  8.  Optimal  point  and  current  curve  as  a  function  of  catalyst  layer 
thickness  L  (Vq  =  0.6  V). 


Fig.  9.  Optimal  point  and  current  curve  as  a  function  of  Nation  volume 
fraction  s n  and  void  volume  fraction  s\  (Vq  =  0.6  V). 
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Fig.  10.  Optimal  point  and  current  curve  as  a  function  of  Nation  volume 
fraction  s n  and  catalyst  layer  thickness  L  (Vq  —  0.6  V). 


Fig.  12.  Optimal  point  and  current  curve  as  a  function  of  Nation  volume 
fraction  s\  and  catalyst  layer  thickness  L  (Vq  =  0.6  V). 


Nafion  and  void  volume  fractions  are  both  taken  into  con¬ 
sideration.  Noting  that  these  two  optimization  parameters 
satisfy  the  constraint  of  Eq.  (17e)  and  this  equation  repre¬ 
sents  a  plane  in  the  space  of  (£n,  £v»  hot)  when  the  volume 
fraction  of  the  solid  catalyst  particles  is  fixed,  therefore  the 
current  density  as  a  function  of  Nafion  volume  fraction  and 
void  volume  fraction  should  be  a  curve  lying  in  the  spatial 
plane  given  by  Eq.  (17e),  as  shown  in  Fig.  9.  The  optimal 
current  density  (0.37467  A/cm2)  occurs  for  a  Nafion  volume 
fraction  and  void  space  volume  of  0.33646  and  0.087743, 
respectively.  These  are  identical  to  the  optimal  values  de¬ 
rived  for  the  individual  one-parameter  cases  shown  in  Figs.  5 
and  6  because  when  one  parameter  is  taken  as  optimizing 
variable,  the  other  varies  according  to  Eq.  (17e)  in  the  opti¬ 
mization  of  each  individual  one-parameter  case. 

Fig.  10  shows  the  current  density  as  a  function  of  the  two 
independent  variables  £n  and  L,  and  the  corresponding  opti¬ 
mal  point  for  the  two-parameter  optimization,  which  occurs 
at  (0.48184,  and  17.5035  ptm)  and  provides  a  maximum  cur¬ 


Fig.  11.  Optimal  point  and  current  curve  as  a  function  of  Nafion  volume 
fraction  s n  and  catalyst  layer  thickness  mpt  (Vq  =  0.6  V). 


rent  0.41186  A/cm2.  Comparing  the  maximum  current  den¬ 
sity  in  Fig.  10  with  those  of  Figs.  5  and  8,  optimization  of 
the  two-parameter  problem  yields  better  results  (higher  cur¬ 
rent)  than  either  of  the  two  one-parameter  optimizations. 

Fig.  1 1  plots  current  density  as  a  function  of  the 
two  independent  variables  £n  and  mpt.  An  optimal 
point  is  observed  from  the  surface  figure  and  is  calcu¬ 
lated  by  the  two-parameter  optimization  as  (0.47627  and 
0.19628  mg/cm2).  Comparing  the  maximum  current  density 
(0.44987  A/cm2)  in  Fig.  11  with  those  in  Figs.  5  and  7,  the 
two-parameter  case  yields  a  much  higher  current  to  the  two 
one- variable  problems. 

The  current  density  as  a  function  of  void  volume  frac¬ 
tion  £y  and  thickness  L  at  a  specified  electrode  potential 
(Vo  =  0.6  V)  is  given  in  Fig.  12.  The  optimal  point  for  (£n, 
L)  is  (0.13927,  and  17.2355  pan)  and  the  maximum  current 
density  is  0.41172  A/cm2. 

The  optimal  values  for  void  volume  fraction  s\  and  Pt 
loading  mpt  are  0.17616  and  0.18641  mg/cm2,  respectively, 


Fig.  13.  Optimal  point  and  current  curve  as  a  function  of  void  volume 
fraction  s\  and  catalyst  layer  thickness  ;npt  (Vq  =  0.6  V). 
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Fig.  14.  Optimal  point  and  current  curve  as  a  function  of  catalyst  layer 
thickness  L  and  catalyst  layer  thickness  (Vo  =  0.6  V). 

and  give  a  maximum  current  density  of  0.45076  A/cm2 
(Fig.  13). 

Fig.  14  gives  the  current  density  as  a  function  of  catalyst 
layer  thickness  L  and  the  platinum  loading  mpt.  A  thickness 
of  2.2807  ptm  and  Pt  loading  of  0.044  mg/cm2  provide  a 
maximum  current  of  0.5 1234  A/cm2  which  is  much  better 
than  those  in  other  cases  above. 

5.  Conclusions 

The  design  optimization  problems  for  cathode  catalyst 
layer  are  discussed  numerically  by  taking  one  or  two  of 
the  four  design  parameters  (sn,  £v>  mvt  and  L)  as  the 
corresponding  optimization  parameter(s).  Numerical  anal¬ 
ysis  shows  that  when  one  of  the  four  design  parameters 
is  taken  as  the  optimization  parameter,  the  corresponding 
one-parameter  optimization  always  has  an  optimal  solu¬ 
tion.  The  catalyst  layer  performance  is  much  more  sen¬ 
sitive  to  the  thickness  than  to  the  other  three  parameters. 
Additional  attention  should  be  paid  to  a  thinner  catalyst 
layer  design  in  order  to  avoid  a  large  drop  in  current 
density.  If  two  of  the  four  design  parameters  are  taken 
together  as  the  optimization  parameters,  the  correspond¬ 


ing  two-parameter  optimization  also  has  a  global  optimal 
solution.  Numerical  analysis  shows  that  by  changing  the 
thickness  and  Pt  loading  the  catalyst  layer  can  reach  its 
best  performance  compared  with  changing  other  pairs  of 
parameters.  The  accuracy  of  the  optimization  problem 
highly  depends  on  the  accuracy  of  the  catalyst  layer  model, 
that  is,  the  mathematical  description  of  the  phenomena 
in  the  catalyst  layer.  Further  theoretical  and  experimental 
studies,  for  example,  the  characterization  of  the  relation 
between  the  effective  and  the  bulk  diffusion  coefficient, 
oxygen  diffusion  in  a  combination  of  Nafion  and  liquid 
water,  should  be  done  in  order  to  improve  the  optimization 
accuracy. 
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